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7—i '. ABSTRACT 

> ; 

' Context. 4U 01 15+63 is one of the most active and best studied Be/X-ray transients. Previous studies of 4U 01 15+63 have led to 

7-H , the suggestion that 4U 01 15+63 undergoes relatively fast quasi-cyclic activity. However, due to the lack of good coverage of the 

■ observations, the variability time scales are uncertain. 

t-H ' Aims. Our objective is to investigate the long-term behaviour of 4U0115+63/V635 Cas to confirm its quasi-cyclic nature and to 

y—i ' explain its correlated optical/IR and X-ray variability. 

, Methods. We have performed optical/IR photometric observations and optical spectroscopic observations of 4U0115+63/V635 Cas 

over the last decade with unprecedented coverage. We have focused on the Ho- line variability and the long-term changes of the 
photometric magnitudes and colours and investigated these changes in correlation with the X-ray activity of the source. 
Results. The optical and infrared emission is characterised by cyclic changes with a period of ~ 5 years. This long-term variability is 
attributed to the state of the circumstellar disc around the Be star companion. Each cycle involves a low state when the disc is very 
weak or absent and the associated low amplitude variability is orbitally modulated and a high state when a perturbed disc precesses, 
5_l ' giving rise to fast and large amplitude photometric changes. X-ray outbursts in 4U 01 15+63 come in pairs, i.e., two in every cycle. 

' However, sometimes the second outburst is missing. 

Conclusions. Our results can be explained within the framework of the decretion disc model. The neutron star acts as the perturbing 
body, truncating and distorting the disc. The first outburst would occur before the disc is strongly perturbed. The second outburst leads 

■ to the dispersal of the disc and marks the end of the perturbed phase. 

Key words, stars: individual: 4U 0115+63, V635 Cas - X-rays: binaries - stars: neutron - stars: binaries close -stars: emission line, 
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1. Introduction 

4U 01 15+63 was one of the first Be/X-ray binaries to be dis- 
covered. The oldest available X-ray observation dates back to 
August 1969 when the Vela 5B satellite dete cted the source 
as thr ee small outbursts separated by 180 days (Whitlock et al. 
119891) . Since then about 15 outbursts have been reported (see 
Table [TJ. Normally, these outbursts represent an increase in 
the X-ray luminosity by a factor ~ 100 and last for about a 
month. The strongest outbursts reach luminosities close to the 
Eddington value (~ 10 38 erg cirT 2 s _1 ). The lack of orbital 
modulation and the relative large increase in luminosity define 
these outbursts as type II. Only in one occasion, in 1996, has 
4U 01 15+63 showed orbitally modulated, s horter and smaller 
outbursts, i.e., type I (Negueruela et al.ll 1998b . For a review of 



* Tables 2 and 3 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/ 



the X-ray variability in Be/X-ray binaries see, e.g., Coe (2000) 
and lZiolkowskH d2002l) . 

X-ray pulsat i ons ( P SD j n = 3.6 s) were soon discovered 
(ICominskv et al.l 1 19781) and the orbital parameters (P 01 -b = 
24.3 d, e = 0.34 , a* sin 2 = 140.1 lt-s) determined 
(Rappaportet al. 19 781). Subsequ ent de tections led to the d iscov- 
ery of one (fWheaton et alii 19791). two dWhite et alJI 19831) . three 
dHeindl et al.ll 19991) and four (ISantangelo et al.ll 19991) cyclotron 
resonance scattering features. 

The 1978 outburst allowed the identification of the opti- 
cal counterpart dJohns et al.ll978tlHutchings & Cramptonl 19811) 
with a reddened (A v > 5 mag) B-type star that showed Ha 
in emission , V635 Cas. Extensive studi e s in t h e IR/optical 
(iKriss et all 119831: iMendelson & Mazehl I199U lUnger etal l 
1 19981: iNegueruela et alj|200ll) and X-rav dWhitlock et al.l ll989: 
iTsunemi & Kitamotolll98 8) bands allowed the determination of 
the astrophysical parameters — V635 Cas is a V ~ 15 B0.2Ve 
star located at a distance of ~ 7-8 kpc — and to the sugges- 
tion of cyclic changes with a quasi-period of 3-5 years. This 
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Table 1. Historical record of the X-ray outbursts and spin period of 4U01 15+63. 



Date 


Satellite 


Intensity 


Reference 


MJD 


p . 

1 spin 


Reference 






(Crab) 










1969 Aug 


Vela 5B 


0.38 


Whitlock et al. (1989) 


40963 


3.614658+0.000036 


Kellev et al. (1981) 


1970 Jan 


Vela 5B 


0.23 


Whitlock et al. (1989) 


42283 


3.6142+0.0001 


Whitlock et al. (1989) 


1970 Aug 


Vela 5B 


0.22 


Whitlock et al. (1989) 


43556 


3.61355+0.000000 


Johns etal. (1978) 


1971 Jan 


Uhuru 


0.17 


Forman et al. (1976) 


44589 


3.6146643+0.0000018 


Ricketts et al. (,198 1 ) 


1974 Aug 


Vela 5B 


1.8 


Whitlock et al. (T989) 


47932 


3.614690+0.000002 


Tamura et al. Q992) 


1978 Jan 


SAS 3 


1.7 


Cominskv etal. (1978) 


47942 


3.61461+0.00001 


Lutovinov et al. (2000) 




Ariel 5 


0.5 


Rose etal. (1979) 


49481 


3.6145107+0.0000010 


Scott etal. (1994) 




HEAO 1 


0.96 


Johns et al. (1978) 


50042 


3.614499+0.000004 


Finger et al. (1995) 


1980 Dec 


Ariel 6 


0.18 


Ricketts et al. (1981) 


50307 


3.614452+0.000003 


Scott etal. (1996) 


1987 Feb 


Ginga 


0.18 


Tsunemi & Kitamoto (1988) 


51232 


3.61452+0.00001 


Wilson et al. (1999) 


1990 Feb 


Ginga 


0.4 


Tamura et al. (1992) 


53261 


3.614+0.003 


Zurita-Heras et al. (2004) 




GRANAT 


0.8 


Lutovinov et al. (2000) 


53266 


3.616+0.001 


Tueller (2004) 


1991 Apr 


CGRO 


0.08 


Cominskv et al. (1994) 








1994 May 


CGRO 


0.18 


Wilson etal. (1994) 








1995 Nov 


CGRO 


0.07 


Finger et al. (1995) 










GRANAT 


0.7 


Sazonov & Sunvaev (1995) 








1999 Mar 


RXTE 


0.4 


This work 










BeppoSAX 


0.3 


Santangelo et al. (1999) 










CGRO 


>0.07 


Wilson etal. (1999) 








2000 Sep 


RXTE 


0.2 


This work 








2004 Sep 


RXTE 


0.3 


This work 










INTEGRAL 


0.4 


Zurita-Heras et al. (2004) 









Vela 5B: 3-12 keV, Uhuru: 2-6 keV, HEAO 1:3-13 keV, SAS 3: 1-27 keV, Ariel 5: 3-6 keV, Ariel 6: 1-50 keV,Ginga: 1-20 keV 
CGRO: 20-50 keV,GRANAT: 8-20 keV, RXTE: 2-12 keV, BeppoSAX: 2-10 keV, INTEGRAL: 20-60 keV 



quasi-cycling behaviour is closely related to the dynamical evo- 
lution of the viscous circumstellar disc around the Be star. The 
Be star loses and reform s the disc on time scales of 3-5 years 
(Negueru ela et alj 12001). At some point during the growing 
phase the disc becomes unstable and highly disturbed. 

In this paper we present the most complete monitoring of 
4U01 15+63/V635 Cas in the optical and infrared wavelength 
bands made up to now. Our observations are used in combination 
with published data to investigate its quasi-cyclic optical/IR/X- 
ray variability. The results presen ted here build on previous work 
by ( Negueruela & Okazaki 2001, Paper I) and lNegueruela et all 
d200Jl Paper II). 

2. Observations 

2.1. Optical and infrared photometry 

Optical (BVRI) photometric observations were obtained from 
the 20cm (before April 2001) and 70cm (after April 2001) tele- 
scopes at the Crimean Observatory (CRI) in Ukraine between 
1998-2005 using an ST-7 CCD of the St.Petersburg University. 
The array is 765 x 5 1 pixels, corresponding to a 8 . 1' x 5 .4' field 
of view (0.65"/pixel scale). The data set consists of about 260 B 
and 350 VRI measurements. In order to improve the signal-to- 
noise ratio, up to 5 images in each colour band were obtained 
and co-added to create the final image. The standard technique 
of bias and dark subtraction and flat-fielding was used. A cali- 
brated set of standard stars, located within the same field, was 
used to perform aperture photometry using a SExtractor-based 
package. The photometric accuracy is usually better than 0^01 
in V, R and / and about 0^03 in B. 

Near-infrared photometric data (about 200 
///^measurements) were obtained at Campo Imperatore 
(Italy) with S WIRCAM NIR camera with the PICNIC 256 x 256 
pixel array, attached to the 1.1m telescope. The camera's field 



of view is 4.5' x 4.5' with a 1.04 "/pixel scale. Each photometric 
image was obtained from 5 co-added dithered images, after 
sky subtraction and flat-field correction. The same photometric 
package as for the ST-7 data was used to perform photometry. 

Table 2 (electronic form only) gives the results of our photo- 
metric monitoring. 

2.2. Optical spectroscopy 

The main source of spectroscopy is the 1.3 m Skinakas telescope 
(SKI), which was equipped with a 2000 x 800 ISA SITe CCD 
and a 1302 1 mrrT 1 grating, giving a nominal dispersion of 1.04 
A/pixel. Many other spectra were obtained through the service 
programme of the Isaac Newton Group at La Palma, either with 
the 2.5 m Isaac Newton Telescope (INT) or the 4.2 m William 
Herschel Telescope (WHT). The INT was equipped with the 
Intermediate Dispersion Spectrograph (IDS), fit with the 235- 
mm camera and different intermediate-resolution gratings. The 
WHT spectra were obtained with the red arm of the Intermediate 
Dispersion Spectroscopic and Imaging System (ISIS), equipped 
with either the R600R or R1200R gratings. 

Other observations have been obtained with the 2.6 m Nordic 
Optical Telescope (NOT), also located in La Palma, equipped 
with ALFOSC; the 1.93 m of the Haute Provence observa- 
tory (OHP) in France, equipped with the Carelec spectrograph 
and the Mt. Ekar 1.82 m Telescope (EKA) of the Padova 
Astronomical Observatory (Italy), equipped with AFOSC. One 
low resolution spectrum was taken with the blue arm of the 
TWIN spectrograph on the 3.5 m of the Calar Alto observatory 
(CA) in Spain. Table 3 (electronic form only) gives the log of the 
optical spectroscopic observations. 

The reductio n of the spectra was ma de using the STARLINK 
Figaro package (Shortridge et al. 2001), while their analysis was 



Reig et al.: Quasicyclic variability in 4U 01 15+63 



3 



1 999 2000 2001 2002 2003 2004 2005 2006 



so 
U 

o 



10 



pi 



> 

PQ 



1 1 

12 

2.5 
2 
1.5 



- T 1 V 1 1 1 ' 

! 

' • »** 

! '» . ••* 

* . 
_ . I ; 

B i tv- »• »*.•*•■" """ 


1 1 1 ~ 

I - 
.. . .^V~* 

*, . «*»•«•• ^- 
* »»«•...> * 

H — " — h- 


H HfH 1 1 1 1 ^ 

K I 

H . r 
J |- •.' fj ' * 

1 ! , 1 


\ •' 

I* *» 

* , 

i 



W 1 

1 

Hi 

0.8 
0.6 



_ 1 

: 

- i 




1 1 ' 1 


1 _ 

^AvU - 


l 1 
- i 


H 1 1 1 1 ' 

[ * 1 , 1 


— 1 — i — H 

1 - 



-20 

"3 

^-10 

UJ 



^800 

1 600 

1400 
< 

200 




1 1 1 1 1 1 1 1 
* 

J 

o 


1 1 1 . 
<f 

CD 

CD 

<t> - 


_ 1 1 1 i 1 J 1 1 1 

o • 

1 j 

o !• < 

> of 4 

_ l ,? 1 , l , l 


; v ; 

1 



51500 52000 52500 53000 53500 54000 

MJD 



Fig. 1. Long-term evolution of the optical and infrared magni- 
tudes and colours and of the Ha equivalent width and peak sepa- 
ration of double-peak profiles. Dotted lines mark the occurrence 
of X-ray outbursts. Filled circles in the bottom panels correspond 
to the observations shown in Fig. [2] 



perfo rmed using the STARLINK Dipso package dHowarfh et alj 
Il998h . 

3. Results 

We have been monitoring 4U 01 15+63 spectroscopically since 
the early 1990's (Paper II). Photometric data were also occasion- 
ally acquired. Since August 1999 the source is being monitored 



in the UBVRIJHK bands. In this section we present the results 
of these new observations. The analysis of the historical variabil- 
ity curves will be presented in Sect. [4] 

3.1. Photometric variability 

Fig.Q]shows the evolution of the optical and infrared magnitudes 
and colours for the period August 2000-June 2006. These obser- 
vations define a complete cycle of variability. The light curve be- 
gins with the photometric bright state that led to the 2000 X-ray 
outburst. The magnitudes gradually decreased and the colours 
became bluer. A faint stable photometric state was reached at 
around MJD 52000 (April 2001). This state extends for more 
than two years up to ~ MJD 52880 (August 2003). We will refer 
to this state as the extended photometric low state (ELS). The 
average BVRI magnitudes during the ELS are B = 16.93 ± 0.05, 
V = 15.50 ± 0.03, R = 14.56 ± 0.03, / = 13.55 + 0.04, 
J = 12.61 + 0.03, H = 12.19 + 0.03, K = 11.86 + 0.05. These 
values are in complete agreement with those given in Paper I for 
the previous faint state. Photometrically, the ELS is a quiet state 
with very low amplitude variations, as can be deduced from the 
small values of the standard deviation. 

The end of the ELS is marked by the gradual brightening 
of the photometric bands in August 2003. The source enters a 
new activity state characterised by large amplitude variations 
in the form of optical eruptions, i.e., the optical/IR brightness 
alternates between maxima and ELS values. The amplitude of 
variability is larger at longer wavelengths (e.g. V increased by 
~0.5 mag and K by 1.6 mag in about 440 days). In coincidence 
with the first optical eruption of this state, a new X-ray outburst 
was observed (September 2004). Note also the constancy of the 
(B - V) colour throughout the observations. 

iBavkal et al.l (120051) reported ROTSE observations covering 
the 2004 outburst. They found a sharp drop lasting for about a 
week (MJD 53235-53242) of ~ 0.3 mag in the source brightness 
a few days before the onset of the X-ray outburst and interpreted 
this result as a sign of mass ejection from the outer parts of the 
disc of the Be star. We do not observe such drop in our data. 
Although our observations contain only one point in the narrow 
interval when the ROTSE drop was seen, it shows very similar 
values to those of the previous and following observations. 

3.2. Spectroscopic variability 

Figure Q] also shows the evolution of the equivalent width 
(EW(Ha')) and the separation of the peaks, in km s~', of the 
split profiles of the Ha line in the interval 1999-2006. As 
the EW(Ha) increases the peak separation decreases. This is 
the expected behaviour of a quasi-Keplerian disk and indi- 
cates that as the EW(Ha') increases, the region where the Ha 
line is produced moves furt her away from the central star (e.g. 
iHummel & Vranckenll995h . 

In Fig. |2] the Ha profile through the different phases of the 
variability cycle is displayed. The top panel corresponds to the 
cycle that began after the 1995 outburst (Paper II), while the bot- 
tom panel to that after the 2000 outburst. The emission line pro- 
file of the Ha line shows significant richness in variability. Both 
single- and double-peak as well as symmetric and asymmetric 
profiles are seen. Symmetric profiles are associated with extreme 
values, either way, of the EW(Ha). Asymmetric profiles are as- 
sociated with intermediate values of EW(Ha). 

Taking as the starting point of the cycle the X-ray outburst 
prior to the ELS, the evolution of the Ha profile through the 
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cycle is the following: at the time of the outburst the Ha line 
shows a single peak profile. The EW(Ha) finds itself in a rela- 
tive maximum. When the profile is a narrow single peak there is 
a broad base component on top of which it sits (see e.g. Fig. [2] 
spectra MJD 51 771 and MJD 51801). T his component was al- 
ready noticed in Neguerue la et al. I (120011) and may be interpreted 
as evidence for a warped disc. The broad base presumably repre- 
sents the motion in the inner part of the disk, where the velocity 
components are higher. After the X-ray outburst the Ha profile 
becomes double peaked and asymmetric with the strength of the 
red peak larger than that of the blue peak (V<R phase). The 
EW(Ha) decreases gradually until the source enters the ELS. 
At this point, EW(Ha) presents its minimum value. Absorption 
as well as double-peak profiles are seen. For some spectra the 
depression between the double-peak profile extends below the 
stellar continuum. As the optical activity increases, the profile 
shows a distorted structure again, although not as marked as dur- 
ing the decay. At the beginning of this phase single-peak and 
red-dominated double-peak profiles are present. As the EW(Ha) 
increases blue-dominated double-peak profiles start to appear. 
When the EW(Ha) reaches ~10 A a new X-ray outburst takes 
place. This outburst does not alter the strength of the Ha line 
significantly, although the overall brightness of the source de- 
creases. By August 2005, the strength of the Ha line and the 
IR emission had experienced a fast and sudden increase. The 
EW(Ha) almost doubled, while the IR magnitudes became ~ 
1 mag brighter, on time scales of a few tens of days. It is as 
though the source was preparing itself for a new X-ray out- 
burst. However, our latest observations (June 2006) show, both 
the EW(Hq') and photometric magnitudes, close to ELS values 
and no evidence of renewed X-ray activity. 



3.3. Reddening and distance 

In Paper I values of the reddening and distance were derived 
from the observations showing the bluest colours (those from 
January 8, 1998). It was then assumed that these observations 
corresponded to purely photospheric emission. In view of the 
observations presented here, a revision of these parameters is 
justified since i) a much longer set of optical photometric obser- 
vations is available, ii) a longer extended low optical state (MJD 
52000-52800) is observed, where presumably the underlying B 
star is exposed and photospheric emission without substantial 
contribution from the disc is detected and Hi) the observations 
cover a wider wavelength band with (quasi)simultaneous optical 
and near IR data, which allows us to estimate the reddening by 
fitting the photometric data to a model a tmosphere. Fig .[3]shows 
the energy distribution of a B0V star dStraizvsl 1 1 9951) and the 
lowest magnitudes of the ELS. Letting Ay be a free parameter 
we find that the best fit is achieved for Ay = 5.17 + 0.03. The 
distributions with Ay = ±3<x are also shown. 

Taking into acco unt that A K = .115Ay, that for a B0V the 
intrinsic Kq - -3.25 (IStraizvsll 19951) and the minimum observed 
magnitude, K - 11.89 ± 0.02, the distance-modulus is DM = 
K - 0.1 15A V - K = 14.55 or d = 8.1 ± 0.1 kpc. The error in the 
distance includes the photometric uncertainty only and assumes 
that neither Kq nor the reddening law are affected by errors — 
for example, an uncertainty of 0.2 mag in Kq translates into a 
distance error of about 0.7 kpc. 

The assumption of standard reddening may introduce large 
errors if not verified. In Paper I, i t was shown that the use of 
the extinction law by Fitzp atrickl d 19991) favours a value of R 
close to the standard R = 3.1. For further verification, we used 



the chorizos code (Maiz-Apellaniz 2004) to study the reddening. 
The p rogram was used to fit extinction laws from lCardelli et all 
dl989b convoluted with the spectral energy distribution of a 
T e s - 27,500 K,logg = 4.0 tlusty atmosphere model to our 
photometry. Although we obtain a slightly lower value for R, the 
derived Ay = 4.8 + 0.1 is compatible with the value found above 
at a 3<x level. 

3.4. Power spectrum analysis 

In order to search for periodic variations in the photometric light 
curves we performed a power spectral analysis. The details of 
the technique employed can be found in lLarionov et al 1 (120011) . 

The power spectrum and spectral window of the data set 
were calculated as 



and 



m)e 



(1) 



w(f) 



N 2 



7=1 



W(0) = 1, 



(2) 



respectively, where i = V— T, tj denotes the time of observation, 
and mj and m, the individual and mean values of the brightness, 
respectively. 

It is clear from Fig. Q] that in order to judge with sufficient 
confidence the reality of the orbital plus close and/or related pe- 
riodicities, one should remove the longer-period component(s) 
of variability. In our case it is not sufficient to subtract the linear 
trend, and we need to take into account the slow light variations. 
We have constructed a smoothed data set using the method of a 
sliding mean with the window value A, replacing raw data m\ for 
each time f; by the weighted mean: 



m i = -2.5 log 



1 



10 



-0.4m j 



(3) 



where k is the number of data points within interval [f; - A, t\ + A], 
and the weight of jth point is determined as 



Pi =exp[-((5fj/A) 2 ], 



(4) 



where 5t s is the time span from the /'th point to the center of 
the window. The optimal value of the smoothing interval was 
searched by trial and error within a range 10 d < A < 30 d . 
The criterion for A selection was the signal-to-noise ratio of the 
power-spectrum peaks in the region of interest, i.e. around the 
orbital frequency. 

We divided the photometric light curves into several intervals 
and analysed them separately: MJD 51400-52000 corresponds 
to the brighter optical state of the source prior to the ELS, MJD 
52100-52850 corresponds to the ELS, MJD 52850-53250 in- 
cludes the smooth brightening after the ELS up to the 2004 X- 
ray outburst and MJD 53250-53715, which covers the time in- 
terval after the X-ray outburst. 

A clear modulation was found during the ELS, independent 
of the value of the smoothing interval adopted. The most promi- 
nent peak in the power spectrum corresponds to a period of 
24.4 + 0.1 d, which is consistent with the orbital period. The 
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Fig. 2. Evolution of the Ha line profile over the ~ 5 year quasi-cycle. The spectra were normalised to the neighbouring contin- 
uum and the wavele ngth converted to velocity units. Indicated are the MJD and orbital phase according to the orbital solution of 
iTamura et all d 19921) . The Y-axis scale was left the same in all panels to facilitate comparison. 



low-frequency components are most effectively suppressed and 
the highest signal-to-noise ratio is obtained when A = 15 d . The 
power spectrum for the residuals Rj - R' is shown in Fig. |4] 

The probability of chance occurrence of a peak with ampli- 
tude P max in the power spectrum with a mean value Pmean can 
be estimated as <D = 100% • {1 - [1 - exp[-(P max /P mean )]] A '">''}; 
for non-uniformly spaced data the number of independent fre- 
quencies Mnd = -6.362 + 1.1 93 ■ N + 0.00098 ■ N 2 , w here N 
is the number of observations dHorne & Baliunaslll986l) . In our 
case N = 301 and Ni n d = 389; after subtraction of the slow com- 
ponent as described above (Eqs.[3]and|4]i, O < 0.2%. 

The amplitude of the sine-wave obtained is 0.01 mag in R. 
We calculate the ephemeris of the small-scale optical variations 
as JDjviax = 2449498(+1.0) + 24.41(+0.15) • E, where E is the 
epoch number. The most striking (but not unexpected) result is 
that these values practically coincide, within errors quoted, with 
the X-ray ephemeris. The position of the optical maximum cor- 
responds to the periastron passage. The orbital period modula- 
tion is also present in the time interval MJD 52850-53250, but 
the shape is markedly non-sinusoidal. 

No evidence of a 24-day peak was found in the interval MJD 
51400-52000. Instead, there is a modulation with 22.15 day pe- 
riod and amplitude « 0.03. One possible explanation for the 
22.15-d period is that it represents the beat frequency between 
the orbital frequency and precession frequency of the disc. In 
this case the disc would precess with retrograde motion with a 
period of about 250 days. After the 2004 X-ray outburst all peri- 
odicities are suppressed. 

We also applied the same analysis to the V and / bands. For 
the / band we obtained practically the same results as in R, while 
in V we were not able to detect any significant modulation close 
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Fig. 3. Energy distribution of a B0V star dereddened using A v = 
5.17. 



to the orbital period. The fact that the beat frequency is more 
apparent in the redder bands gives support to the interpretation 
of the 22.15-d periodicity as being affected by the precession 
of the disc, as one would expect the IR magnitudes to be more 
affected by the disc that the shorter wavelengths. 
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Fig. 4. Top: R-R' dependence on the phase of the 24. 31 -day pe- 
riod for the ELS. Bottom: Power spectrum of the residuals R-R' 
obtained from the period MJD 52100-52300. The most promi- 
nent peak corresponds to P = 24.41 days. 



4. Discussion 

4.1. Source states 

We have analysed optical light curves and spectra of 
4U01 15+63 spanning over more than two decades, the last 
10 years with unprecedented coverage. During this period the 
source exhibited significant photometric and spectroscopic vari- 
ability, allowing the definition of source states. The intensity of 
the optical, infrared and X-ray emission highly correlates with 
the state of the source. 

Figure [5] shows the historical record (1980-2006) of the R 
band and Ha equivalent width (EW(Ha)) of V635 Cas. Squares 
and stars i n the top panel of Fig. [5] repre sent observations 
taken from iMendelson & Mazehl dl99ll) and iNegueruela et al.l 
( 2001b. respectiv e ly. Sq uares in the EW( Ha) panel a r e from 
INegueruela et al.l d2001|) and crosses from iKriss et al.l {[983). 
Circles correspond to our new measurements. Vertical dashed 
lines denote the occurrence of X-ray outbursts. Note the strik- 
ing repeatability of the pattern of variability during the periods 
1986-1990 and 2002-2005: a photometric faint sate is followed 
by an optical flare with a slow gradual rise (200-300 days) and 
a fast decay (70-80 days). An X-ray outburst occurs when the 
source is in a photometric bright state. A second optical eruption 
separated from the previous one by about 250 days is not accom- 
panied by X-ray activity. This second optical eruption presents 
a m ore symmetric profile with a fast rise and equally fast de- 
cay. |Mendekon]&^Mazeh| {l99l]) also reported the presence of a 
small flare after the second eruption, peaking in June 1988 (MJD 
47340) with an amplitude of ~ 0.3 mag and duration of 50 days. 
Exactly the same event is seen 17 years after (MJD 53585). If 
the source follows the same pattern as in the late 1980's, then 
we should expect the source to remain in a low photometric state 
until spring 2007 and to start a new optical outburst soon after. 
The overall difference of ~ 0.25 mag between the Mendelson 
& Mazeh's and our data set is most likely due to instrumental 
as well as calibration (selection of different secondary standard 
stars) e ffects, as no attempt to perform absolute photometry was 
made in IMendelson & M azeh (1991). The optical/IR/X-ray be- 
haviour of 4U0115+63/V635 Cas can be understood in terms 



of the evolution of the equatorial disc around the Be star. Two 
basic states can be distinguished depending on the presence or 
absence or the disc. 



AAA. The low state 

The low state would correspond to the complete loss of the 
disc or to a highly debilitated disc. During the low state, the 
source shows the weakest magnitudes, the bluest colours and 
the smallest Ha equivalent widths (EW(Ha)). The source ex- 
hibits little variability, with changes in the photometric bands of 
less than 0.05 magnitudes. Orbital modulation is detected. Since 
the source stays in the low state for extended peri ods (a few 
years ), the term 'extended low state' (ELS) is used dRoche et al.l 
1993). During the low state the source is X-ray quiet. A disc- 
less state occurs when the Ha line displays an absorption pro- 
file. According to our interpretation of the long-term light curve 
shown in Fig. [6] 4U 01 15+63 would have gone through ELS dur- 
ing MJD 48600-49300 (1991-1993), MJD 50300-50900(1996- 
1998) and MJD 52000-52800(2001-2003). 

The issue of whether the disc completely vanishes during 
the low state or there still exists some residual emission from a 
highly debilitated disc is a very important one, since in the case 
of total loss one can decouple the disc emission from that of the 
central star and hence determine the astrophysical parameters of 
the underlying star without any interference from the disc (Paper 
I). Disc-less states in 4U01 15+63 occurred in 1997 and 2001 as 
the Ha line appeared in absorption. Further support for the com- 
plete disappearance of the disc in 2001 comes from the fact that 
the standard deviation of the photometric magnitudes during this 
state does not follow any trend as a function of wavelength (see 
Sect. 13.11 ). Disc activity is expected to affect the redder magni- 
tudes more than bluer magnitudes, hence making the IR magni- 
tudes appear more variable. Also, from the power spectral anal- 
ysis we observed that when the disc is present a beat period be- 
tween the orbital period and presumably the precession period of 
the disc was found. Then the fact that no such beat frequency is 
detected during the ELS but only a low-amplitude (~ 0.01 mag) 
modulation coinciding with the orbital period seems to indicate 
the absence of a disc. The 1992 low state would have not been 
accompanied by the complete loss of the disc as no absorption 
profile was attained. 

The neutron star spins down during the low state. Previous 
studies have reported spin-up e pisodes throughout t he duration 
of the X-ray outbursts (see e.g iTamura et alj |1992h . However, 
the historical record of the spin period hardly shows evidence 
for variability (see Table [1]). Thus, if the neutron star spins up 
during the outburst, it must spin down in quiescence. 

Disc-less phases have been seen in a number of Be/X-ray 
binaries but only two have been observed frequently enough as 
to allow a me aningful comparison wit h 4U 01 15+63, namely , 
A 0535+26 dLvutv & Zaitseval 120001; lLarionov et all [2 001: 
[Clark et all Tl 9991: lHaigh et all 120041) and X Per dClarketal.1 
1200 ll and references therein). The three systems show extended 
low states, although its duration differs: -600 days (MJD 5 0700- 
51300) for A 0535+26 dHaigh et al.ll200l:IZaitsevall2005l) . -800 
days (MJD 52000-52800) for 4U0115 +63 (Fig. □ > and - 1300 
days (MJD 47700^9000) for X Per dClark et alj|200lb . The 



long-term photometric light curv e of X Pe r is ve ry similar to 
that of 4U 01 15+63 (cf. Fig 7 in lClark et al.l d200l with FigQ), 
whereas A 0535+26 ex hibits a more choppy behaviour (see 
Fig. 1 in lClark et alJfl999l) . 
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Fig. 5. Long-term evolution o f the R band and the EW(Hcr ). In 
the top pa nel squares are fromlMendelson & Mazehl dl99lh and 
stars from Neguerue la et al.l (2001). In the bottom pa nel squares 
are fro m [Neguerue la et al.l (12001) and crosses from Rris s et all 
d!983l) . The data from the year 2000 onward represent new ob- 
servations. Dotted lines mark the occurrence of X-ray outbursts. 



4.1.2. The high state 

As the disc reforms, it first behaves like the disc of an isolated 
Be star, growing in size. This phase is relatively smooth and it is 
characterised by a gradual increase in brightness and EW(Ha). 
After the disc has reached its maximum size, further mass loss 
will result in an increase of the disc density, which can even 
become optically thick at IR wavelengths. At this point, the con- 
tribution of the disc to the optical/IR photometric magnitudes 
and colours is substantial. Fast (~100 days) and large amplitude 
variations (upto 1 .5 magnitudes) are seen. The source brightness 
alternates between low-state values and optical maxima. In this 
state the shape of the Ha line alternates between single-peak and 
double-peak profiles, indicative of a perturbed disc. 

X-ray outbursts always occur during the high state. They are 
preceded by the brightening of the optical and infrared emission. 
The EW(Ha') reaches a relative maximum at the time of the X- 
ray peak. There seems to exist a triggering value of the EW(Ha), 
around 10 A, for the onset of the X-ray outburst. The EW(Har) 
prior to the start of the second outburst is always larger than 
that of the first one. Each X-ray outburst is then followed by a 
decrease in the optical/IR brightness but not always by the loss 
of the disc. Note that the EW(Ho') remained fairly unaltered after 
the 1994, 1998 and 2004 X-ray outbursts (Fig.©. 

When the X-ray outburst does not lead to a low state (i.e. disc 
loss) the EW(Ho') presents a relative minimum a few months af- 
ter the outburst (Table 3, see also Table 1 in Paper II), as if the 
disc would have started a dispersal phase that never materialises. 
If this disc-loss phase does not take place, then another X-ray 
outburst occurs about 1—1.5 years later and is promptly followed 
by the complete disappearance of the disc. Indeed, a very inter- 
esting feature of the long-term light curve of 4U01 15+63 is the 
fact that, for the last two decades, X-ray outbursts appear to have 
occurred in pairs, with two major outbursts happening in every 
cycle of disc build-up and loss (Fig.©. 

Although a comparison of the strength of the outbursts is 
hampered by the different energy ranges covered by the instru- 
ments that detected the events, we observe that the X-ray lumi- 
nosity of the first outburst tends to be larger than the second one 
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Fig. 6. Our interpretation of the long-term light curve of V635 
Cas. Marked points represent real observations. Low states occur 
in between pairs of X-ray outbursts, high states after the first 
outburst. 



(Table©. In one occasion (in 1996, MJD 50300-50400) the sec- 
ond X-ray outburst was followed b y a series of orbit al modulated 
minor outburst (so-called type I) (Negueruela et al. 1998). 

The second outburst of the pair, however, is not always de- 
tected. This seems to be the case in the last cycle (2000-2005). 
Photometric magnitudes are back to low state values and 1.5 
years have passed since the September 2004 outburst. Our latest 

spectra taken on June 2006 give values of the EW(Hq') of 0.4 

A and an absorption profile with two weak emission shoulders, 
a feature typically seen in small disks. This, however, does not 
necesarily mean that 4U01 15+63 has returned to the low state. 
In Paper II we showed that fast transitions between single-peak 
and shell profiles are indicative of a precessing warped disc. 
Another missing outburst oc curred during the 1985-1990 cycle 
dMendelson & Mazehlll99ll) . The similarity between the 2003- 
2006 and 1986-1989 light curves is striking. The second optical 
eruption in 1988 was not accompanied by an X-ray outburst ei- 
ther. Therefore we suggest that the intervals MJD 46700-48600 
(1986-1991) and MJD 52800-present (2003-present) should be 
considered as extended high states. Extended high states would 
occur in association with missing X-ray outbursts. Normal high 
states would last for about 3 years — the intervals MJD 47600- 
48600 (1989-1991), MJD 49300-50300 (1994-1996) and MJD 
50900-52000 (1998-2000) would correspond to high states — 
while extended ones would cover the entire cycle of 5 years. In 
this case no disc-less phase is achieved. 

4.2. The viscous decretion disc model 

The long-term variability of 4U 01 15+63 can be summarised as 
follows: the equatorial disc around the Be star dictates the state 
of the source. The quasi-cyclic variations coincide with the time 
scales for disc build-up and loss. X-ray outbursts occur in pairs, 
although the second one may be missing. In between pairs of 
outbursts the optical/IR emission and variability is highly sup- 
pressed, while in between the outbursts of a given pair the opti- 
cal/IR emission is highly variable. 

This behaviour can be explained wit hin the framework of the 
trunc a ted decretion disc m odel (Paper I, Okazaki & Neg ueruelal 
1200 U lOkazakiet ail 120021) . Initially, the disc grows gradually 
like in isolated Be stars. However, in Be/X-ray binaries the 
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surface density of the disc increases more rapidly that that of 
isolated Be stars, as a consequence of truncation. The trun- 
cation of the disc is the result of the resonant torque exerted 
by the n eutron star, which re moves angular momentum from 
the disc. lOkazaki et ail d2002l) locates the truncation radius of 
4U 01 15+63 at 0.36 - 0.39a, where a is the semi-major axis 
of the orbit, for a wide range of values of the viscosity. The 
disc becomes optically thick at IR wavelengths and unstable to 
radiation-driven warping. Eventually, the disc begins to warp, tilt 
and precess. 

The photometric variations observed during the high states 
(1987-1989, 2003-2005), characterised by large (> 1 mag) and 
relatively fast (~100 days) amplitude variations, are rather large 
when compared t o typical values seen in isolated Be stars 
dDachs et al.ll9 88). In addition, the brightest magnitudes are ac- 
companied by historic maxima of the EW(Ha). These extreme 
values of the EW(Ha) correspond to line shapes much narrower 
than at other times. In the framework of the disc deccretion 
model, the coincidence of very high values of the EW(Ha) with 
narrow, single-peaked lines and very bright optical and infrared 
magnitudes during the high state, leads to the conclusion that 
extreme values happen because the disk becomes warped and 
occasionally presents a much higher surface to the observer (see 
Paper II). Therefore the changes in magnitudes should be the 
consequence of two different effects: on the one hand, an in- 
crease in the disk size or density; on the other hand, a change in 
the emitting surface observed. They should not be interpreted as 
fast phases of disc loss and reformation, despite the fact that the 
photometric magnitudes fall to ELS values. This strongly dis- 
turbed phase would begin after the first outburst. 

One possible explanation for the difference between the first 
and second outburst could be the different state of the disc, with 
the first one occurring when the disc is still in a quasi-stable state 
while the second one when the disc has been distorted. The dis- 
tortion of the disc leads to the interaction with the orbiting neu- 
tron star. If the interpretation of the 22.15 day period (Sect. [3~4l 
as the beat period between the orbital period and precession is 
correct, then the time scales associated with the perturbed disc 
(~250 d) would be a factor ~2 longer than the duration of the op- 
tical eruption associated with the missing X-ray outburst (~150 
d). In this context the missing outburst would simply reflect the 
fact that no interaction between the neutron star and the distorted 
disc took place during the ~6 orbits that the optical eruption 
lasted. 

An alternative explanation for the missing X-ray outburst 
requi res the formation of an acc retion disc around the neutron 
star dMend elson & Mazeh ¥T99 ll) and it is based on the centrifu- 
gal inhibition at the edge of the neutron star magnetosphere 
dStellaet al.lll986l) . In the 1988 and 2005 events the matter ac- 
cumulated in the disc would not be enough to overcome the pro- 
peller effect. 

The (B-V) colour index appears as a good indicator to distin- 
guish between the two outbursts. While the (J - K) colour (also 
(V—R) and (R-I), not shown in Fig[T|i follows smoothly the state 
of the disc, (B - V) is rather insensitive to photometric or spec- 
tral changes. Only prior to the destructive outburst of September 
2000 does it display high amplitude variability (Fig. |T}. 

In addition to the missing outbursts, the fact that the first X- 
ray outburst does not modify, at least dramatically, the strength 
of the Ha line also needs to be explained. In this respect, one 
may wonder how much of a disruption an X-ray outburst repre- 



sents. As an order of magnitude estimate, we can derive the mass 
content in the disc simply as 

M = j pdV 

with p = po{R*/r)" and dV = InrHdr. Here po is the density 
and the inner radius of the disc (at r ~ R„), R, the star radius, 
H the disc hei ght and n an expon ent defining the density law. 
According to Wate rs et al.l JT988) n varies in the range 2-4 in 
most BeX. The typical radius of a B0. 2 star, is 7?« = 8 R and 
assuming typical val ues for the disc rad ius R& = 5R t , and inner 
density po = 10" 10 g (Telling etalJl 1998b cm 3 and H = 0.03R* 
dNegueruela & Okaz aki 2001) the mass is estimated to be 2.6 x 
10~ 9 M Q for n = 2 and 7.8 x 10~ 10 M Q for n = 4. 

On the other hand, a mean X-ray luminosity of 10 37 erg s _1 
for about a month requires the transfer of M — L X {R/GM) = 
5.4 x 10 16 g s" 1 or 7 x 10~ n M month -1 , that is just < 10% 
of the mass in the disc. This relatively small effect on the mass 
content of the disc would explain the fact that the strength of the 
Ha line does not c hange much after the first outburst (see also 
iNorton et alJll994h . 

5. Conclusion 

We have presented the results of our monitoring of the Be/X- 
ray binary 4U0115+63/V635 Cas. Our 2000-2006 data repre- 
sent the most complete optical and infrared photometric study 
of 4U 01 15+63/V635 Cas made up to now. The combination of 
our observations over the last decade with published data al- 
lowed us to investigate the correlation between the optical, in- 
frared and X-ray emission. X-ray outbursts, which come into 
pairs, occur when the source is photometrically bright and the 
Ha line appears strongly in emission. The equatorial disc around 
the Be star dictates the state of the source: if the disc is absent 
the source is in the low state. In the high state large amplitude 
variation and asymmetric spectral lines denotes a disturbed disc. 
The time scale for loss and build-up of the equatorial disc is 
5 years, approximately 3 of which correspond to the high state 
and 2 to the low state. However, this quasi-periodic behaviour is 
broken if one of the X-ray outburst is missing. In this situation 
the high state covers the entire cycle. Bright magnitudes are ob- 
served because the disc warps and presents a larger surface to 
the observer. The fact that the infrared and optical magnitudes 
and the Ha equivalent width do not change dramatically after 
the first outburst may indicate that, at the time of the outburst, 
the disc had not been distorted yet. It is after the first outburst 
that the disc becomes unstable, warps and tilts. The second X- 
ray outburst takes place during this phase of strongly disturbed 
disc. The disc finally disappears, partly reabsorbed by the B star 
and partly used up to power the X-ray outburst. 

Acknowledgements. We thank A. Manousakis, A. Di Paola and M. Dolci for 
their help in some of the observations. IN is a researcher of the programme 
Ramon y Cajal, funded by the Spanish Ministerio de Educacion y Ciencia 
(MEC) and the University of Alicante, with partial support from the Generalitat 
Valenciana and the European Regional Development Fund (ERDF/FEDER). 
This research is partially supported by the MEC under grant AYA2005-00095. 
VL and AA acknowledge hospitality of the Rome and Teramo observatories dur- 
ing the observations in Campo Imperatore. This work has greatly benefited from 
the ING service programme. Skinakas Observatory is a collaborative project of 
the University of Crete, the Foundation for Research and Technology-Hellas 
and the Max-Planck-Institut fur Extraterrestrische Physik. The near-infrared 
data were obtained with AZT-24 telescope operated under agreement between 
Pulkovo, Rome and Teramo observatories. The INT and IKT are operated on 
the island of La Palma by the Isaac Newton Group in the Spanish Observatorio 
del Roque de Los Muchachos of the Instituto de Astroffsica de Canarias. This 



Reig et al.: Quasicyclic variability in 4U 01 15+63 



9 



research is partially based on data from the ING Archive. Part of the data pre- 
sented here have been taken using ALFOSC, which is owned by the Institute de 
Astrofisica de Andaluci'a (IAA) and operated at the Nordic Optical Telescope un- 
der agreement between IAA and the NBIfAFG of the Astronomical Observatory 
of Copenhagen. Partially based on observations obtained at the Observatoire de 
Haute Provence (CNRS, France), observations obtained at the Asiago obser- 
vatory and the Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, 
operated jointly by the Max-Planck Institut fr Astronomie and the Instituto 
de Astrofisica de Andaluci'a (CSIC). This research has made use of NASA's 
Astrophysics Data System Bibliographic Services and of the SIMBAD database, 
operated at the CDS, Strasbourg, France. 



References 

Baykal, A., Kiziloglu, U., Kiziloglu, N., Balman, S., & Inam, S. C. 2005, A&A, 
439,1131 

Cardelli, J.A., Clayton, G.C., Mathis, J.S. 1989, ApJ, 345, 245 

Clark, J. S., Lyuty, V. M., Zaitseva, G.V., et al., 1999, MNRAS, 302, 167 

Clark, J. S., Tarasov, A. E., Okazaki, A. T, Roche, P., Lyuty, V. M. 2001, A&A, 

380, 615 
Coe, MJ 2000, ASPC, 214, 656 

Cominsky, L., Clark, G. W., Li, F., Mayer, W., & Rappaport, S. 1978, Nature, 
273, 367 

Cominsky, L., Roberts, M., & Finger, M.H. 1994, in The Second Compton 
Symposium, AIP Conf. Proc. 304, Ed C. E. Fichtel, N. Gehrels, J. P. Norris, 
AIP Press, New York, NY, p. 294-298. 

Dachs, J., Engels, D., Kiehling, R. 1988, A&A, 194, 167 

Finger, M. H., Scott, M., Hagedon, K., Wilson, R. B., Wilson, C. A., Stollberg, 

M., Paciesas, W. S., Prince, T. A., Vaughan, B. 1995, IAUC, 6266 
Fitzpatrick, E.L., 1999, PASP 111, 63 

Forman, W., Jones, C, & Tananbaum, H. 1976, ApJ, 208, 849 

Haigh N.J., Coe M.J., Fabregat J. 2004, MNRAS, 350, 1457 

Heindl, W. A., Coburn, W., Gruber, D. E., Pelling, M. R., Rothschild, R. E., 

Wilms, J., Pottschmidt, K., & Staubert, R. 1999, ApJ, 521, L49 
Home J.H., Baliunas S.L., 1986, ApJ 302, 757 

Howarth, I. D., Murray, J., Mills, D., & Berry, D. S. 1998, Starlink User Note 
50.21 

Hummel, W. & Vrancken, M. 1995, A&A, 302,751 
Hutchings, J. B. & Crampton, D. 1981, ApJ, 247, 222 
Johns, M„ Koski, A., Canizares, C, & McClintock, J. 1978, IAUC 3171 
Kelley, R. L., Brodheim, M. J., Cominsky, L., Stothers, R., Rappaport, S. 1981, 
ApJ, 251, 630 

Kriss, G. A., Cominsky, L. R., Remillard, R. A., Williams, G, & Thorstensen, J. 

R. 1983, ApJ, 266, 806 
Larionov, V., Lyuty, V. M., & Zaitseva, G. V. 2001, A&A, 378, 837 
Lutovinov, A. A., Grebenev, S. A., & Sunyaev, R. A. 2000, AstL, 26, 1 
Lyuty V.M., Zaitseva GV. 2000, AstL, 26, 9 
Maiz-Apellaniz, J. 2004, PASP, 116, 859 
Mendelson H, Mazeh T. 1991, MNRAS, 250, 373 
Negueruela, I., Reig, P., Coe, M. J., Fabregat, J., 1998, A&A, 336, 251 
Negueruela, I. & Okazaki, A. T. 2001, A&A, 369, 108, Paper I 
Negueruela, I., Okazaki, A. T, Fabregat, J., Coe, M. J., Munari, U., Tomov, T. 

2001, A&A, 369, 1 17, Paper II 
Norton, A. J., Chakrabarty, D., Coe, M. J. et al. 1994, MNRAS, 271, 981 
Okazaki, A.,T. & Negueruela, I. 2001, A&A, 377, 161 

Okazaki, A.T., Bate, M.R., Olgivie, G.I., & Pringle, J.E. 2002, MNRAS, 337, 
967 

Rappaport, S., Clark, G. W., Cominsky, L., Li, R, & Joss, P. C. 1978, ApJ, 224, 
LI 

Ricketts, M. J., Hall, R., Page, C. G, & Pounds, K. A. 1981, SSRv, 30, 399 

Roche, P., Coe, M. J., Fabregat, J., et al. 1993, A&A, 270, 122 

Rose, L. A., Marshall, F. E., Holt, S. S., Boldt, E. A., Rothschild, R. E., 

Serlemitsos, P. J., Pravdo, S. H, Kaluzienski, L. J. 1979, ApJ, 231, 919 
Santangelo, A., Segreto, A., & Giarrusso, S. 1999, ApJ, 523, L85 
Sazonov S., & Sunyaev R. 1995, IAUC. 6272 

Scott, D. M., Finger, M. H, Wilson, R. B., Prince, T. A. 1994, IAUC 5990 
Scott, M., Finger, M. H, Wilson, R. B., Prince, T. A., Vaughan, B. 1996, IAUC 
6450 

Shortridge, K., Meyerdierks, H, Currie, M., et al. 2001, Starlink User Note 86.19 
Stella, L., White, N. E., Rosner, R. 1986, ApJ, 308, 669 

Straizys V., 1995, Multicolor stellar photometry, Astronomy and Astrophysics 
Series, vol.15, Parchart Publishing House, Tucson, Arizona 

Tamura, K., Tsunemi, H., Kitamoto, S., Hayashida, K., Nagase, F. 1992, ApJ, 
389,676 

Telting, J.H., Waters, L.B.F.M., Roche, P., Boogert, A.C.A., Clark, J.S., de 

Martino, D., Persi, P. 1998, A&A, 339, 150 
Tsunemi, K. & Kitamoto, S. 1988, ApJ, 334, L21 



Tueller, J. 2004, IAUC 841 1 

Waters, L.B.F.M., van den Heuvel, E.P.J., Taylor, A.R., Habets, G.M.H.J., Persi, 

P. 1988, A&A, 198, 200 
Wheaton, W. A., Doty, J. P., Primini, F. A., Cooke, B. A., Dobson, C. A., 

Goldman, A., Hecht, M., Howe, S. K., Hoffman, J. A., Scheepmaker, A. 1979, 

Natur, 282, 240 
White N.E., Swank J.H., Holt S.S 1983, ApJ, 270, 711 
Whitlock, L., Roussel-Dupre, D., & Priedhorsky, W. 1989, ApJ, 338, 381 
Wilson, R. B., Finger, M. H, & Scott, D. M. 1994, IAUC 5999 
Wilson R.B., Harmon B.A., Finger M.H. 1999, ATel 40 

Unger, S. J., Roche, P., Negueruela, I., Ringwald, F. A., Lloyd, C, Coe, M. J. 

1998, A&A, 336, 960 
Zaitseva GV. 2005, AstL, 31, 103 
Ziolkowski, J. 2002, MmSAI, 73, 1038 

Zurita Heras, J. A., Shaw, S. E., Kreykenbohm, I. et al. 2004, ATel, 331 



